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Abstract 

We propose a consistent scenario of the evolution of the universe based on 
the large cutoff supergravity (LCSUGRA) hypothesis of supersymmetry breaking, 
where the gravitino and sfermion become as heavy as ~ 0(1 — 10 TeV). With such 
a heavy gravitino, baryon asymmetry of the universe can be generated by the non- 
thermal leptogenesis via an inflaton decay without conflicting the serious gravitino 
problem. We also see that, in the LCSUGRA scenario, relic density of the lightest 
superparticle becomes consistent with the WMAP value of the dark matter density 
in the parameter region required for the successful non-thermal leptogenesis. 



1 Introduction 



In a recent paper pQ Izawa and two of us (M.I. and T.Y.) have proposed a large-cutoff 
hypothesis in supergravity. In this hypothesis all higher dimensional operators such as 
quartic terms in the Kahler potential at the GUT scale Mqut — 2 x 10 16 GeV (or at the re- 
duced Planck scale M G ~ 2.4 x 10 18 GeV) are suppressed by a large cutoff M* ~ 0{AnM G ). 
Then, the sfermions and gravitino become order-of-magnitude heavier than the gauginos 
and, consequently, masses of the sfermions and gravitino are required to be significantly 
larger than the electroweak scale. Even so, naturalness of the electroweak symmetry 
breaking can be maintained by the focus-point mechanism [2] due to the fact that the 
universality of the scalar masses at the GUT scale is guaranteed in this scenario. In pQ, 
it was shown that the large cutoff supergravity (LCSUGRA) scenario is well consistent 
with low-energy phenomenology. In particular, heaviness and universality of the sfermion 
masses are good for suppressing dangerous supersymmetric effects on the flavor violating 
processes, proton decay, and so on. We consider that the presence of the large cutoff M* 
is a reflection of a more fundamental physics beyond the GUT scale. 

In this letter we study the cosmology of the LCSUGRA scenario. In SUGRA, there is a 
serious cosmological problem, that is the gravitino problem. If the gravitino is unstable, it 
has a long lifetime and decays after the big-bang nucleosynthesis (BBN) for an interesting 
range of the gravitino mass, 7713/2 ~ lOOGeV — lOTeV. The decay products destroy light 
elements produced by the BBN and hence the primordial abundance of the gravitino is 
constrained from above to keep the success of the BBN. This leads to an upper bound 
on the reheating temperature Tr after inflation, since the abundance of the gravitino 
is proportional to Tr. The recent detailed analysis derived a stringent upper bound 
Tr^S 10 6 ~ 7 GeV when the gravitino decay has hadronic modes jSj- This upper bound is 
much lower than the temperatures required for most of thermal baryogenesis including 
the leptogenesis [I]. There have been proposed various solutions to the above problem 
and, among them, we consider that the non-thermal leptogenesis via an inflaton decay 
jS] is the most interesting and plausible. 1 Importantly, even if the right-handed neutrinos 

1 The Affleck-Dine (AD) baryogenesis (j is an interesting scenario, but this does not work in the 
LCSUGRA, since the Kahler potential between inflatons and the AD fields is suppressed and the AD 
fields acquire the Hubble- induced masses which set the AD fields at the origin during the inflation jxj. 
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are non-thermally produced from the decay of the inflaton, success of the baryogenesis is 
not automatic, as we will see. 

The main conclusion of this letter is that the LCSUGRA scenario is also advantageous 
in cosmology In the next section, we show that a general argument on the inflaton- 
decay scenario of leptogenesis gives us a lower bound on the reheating temperature as 
Tr^jI x 10 6 GeV jS]. Together with the constraint from the BBN [3] we find the lower 
bound on the gravitino mass m.3/2 ^ 4 TeV. Importantly, such heavy gravitino is a natural 
outcome of the LCSUGRA scenario. In sectional we show that the LCSUGRA with such 
a large gravitino mass indeed explains the observed dark matter density. In particular, 
the result of the present analysis shows that the universal scalar mass of SUSY breaking, 
mo, or equivalently the gravitino mass 711,3/2 is 4 — 10 TeV which is in a special area of 
so-called focus point region. Remarkably, such a parameter region is also required for the 
success of the non-thermal leptogenesis via the inflaton decay as mentioned above. The 
last section is devoted to discussion. 

2 Leptogenesis via inflaton decay 

Let us first consider the leptogenesis in the framework of LCSUGRA. Here, for simplicity, 
we consider the case where only the lightest right-handed neutrino N\ contributes to the 
lepton asymmetry, assuming mass hierarchy among the right-handed neutrinos. N% may 
decay into H u + £ or H* + £* where H u is the up-type Higgs and £ the lepton doublet. 
These two decay channels have different branching ratios when CP conservation is vi- 
olated. Interference between tree-level and one-loop diagrams generates lepton-number 
asymmetry jU HU] as 

r(N 1 ^H u + e)-r(N 1 ^H* u + t) 3 m 1 

T Nl 8% {H u } 2 

where m Vz is the heaviest (active) neutrino mass. Using 3x3 neutrino Yukawa matrix h, 
the effective CP-violating phase <5 e fr is given by 
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In deriving above expressions, the seesaw mass formula ^1] has been used. 



3 



With non-vanishing e parameter, lepton-number asymmetry is generated with the 
out-of-equilibrium decay of N^. Since e is proportional to M±, Nx is required to be 
heavy enough to generate sufficient amount of baryon number density. If we consider 
the case where the right-handed neutrino is thermally produced, for example, we obtain 
Mi> 10 9 - 10 GeV im Since Tr ^ M\ in this case, we should conclude that the thermal 
leptogenesis requires too high reheating temperature to avoid the gravitino problem in 
the conventional SUGRA models. 

If the right-handed neutrino is non-thermally produced by the decay of the inflaton 
$ as $ — > N\ + N\ ($ — > Ni + Ni) [5], the situation changes. In this case, Boltzmann 
equations for the B — L number density ub-l and other quantities (i.e., the energy density 
of the inflaton p$ and that of radiation p ra d) are given by 



driB-L 
dt 

dp<$, 
~dt 

dprad 



+ 3Hn B - L = 2er$m # 1 p* ) ( 3 ) 
3Hp$ = -r$p$, (4) 



+ AHp rad = T^ps, (5) 

where m$ is the inflaton mass, T$ the decay rate of the inflaton. Notice that the factor 
2 in the right-hand side of Eq. (JSJ) is due to the fact that two right-handed neutrinos are 
produced by the decay of single Here and hereafter, we adopt a mild assumption that 
the inflaton potential can be well approximated by the parabolic one at the last stage of 
the reheating. 

The inflaton decays when the expansion rate of the universe becomes comparable to 
the decay rate we define the reheating temperature as 

where g* is the effective number of the massless degrees of freedom. (In our numerical 
study, we use g* = 228.75.) As well as the entropy production, generation of the B — L 
asymmetry is most effective when H ~ Consequently, we obtain the relation between 
ub-l and the entropy density s as 

n B-L T R n / T R \ (2M X \ ( m U3 \ 

= Ke = 9.9 x 10 x ft — — — — <W, (7) 

s m $ V10 6 GeV/ V m $ A0.05eV/ ' y J 
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where k is a constant of 0(1). In the second equality, we have taken (H u ) ~ 174 GeV. 
Qualitative behavior of Eq. (JJJ) can be easily understood using the instantaneous-decay 
approximation: m^riB-LiTR) ~ e[p^]deca y ~ e 9*T^ and s(Tr) ~ g*T%- Using the order-of- 
magnitude relation tib-l ~ after the spharelon transition, one can see that, in order 
to generate the enough baryon asymmetry of the universe suggested by the WMAP [T3] 

^~0.9xl(T 10 , (8) 

Tr should be higher than ~ 10 6 GeV. 

Since this lower bound is very close to the upper bound on Tr from the gravitino 
problem, we have performed a careful calculation of the resultant baryon number in this 
scenario. We have numerically solved Eqs. © — © from the cosmic time t <C Tj 1 to 
t ^> Tj 1 . Then, we have calculated the entropy density at t ^> Tj 1 using the relations 

Taking the ratio of ub-l to s, we have obtained the k parameter in Eq. © as 

k ~ 2.44. (10) 
Using ub = jk n B~L [T4l . baryon-to-entropy ratio is given by 



n B 
s 



8.2 x 10-" x (**-) (™L) (J**-) 5 eS . (11) 
U0 6 GeV A m$ Ao.05eW V ; 

Combining this with Eq. ©, we derive a constraint on the reheating temperature, 

T R >lxl0 6 GeV, (12) 

for the neutrino mass suggested from the atmospheric neutrino oscillation experiments, 
Ara| 3 — 0-05 eV. Here, we used Mi = |m$, which is the maximal possible value 



of Mi so that the decay process $ — > N± + N\ is kinematically allowed. From the point 
of view of the gravitino problem, such reheating temperature is quite dangerous. 

It is notable that, in the LCSUGRA scenario where the gravitino is quite heavy, the 
present non-thermal leptogenesis scenario becomes viable in a wide parameter region. 
Indeed, the recent detailed analysis of the gravitino problem suggests that Tr ~ 1 x 
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10 6 GeV is consistent with the BBN constraints if the gravitino is heavier than 4 TeV [Sj. 
In the LCSUGRA scenario, such a large gravitino mass can be realized without conflicting 
the naturalness of the electroweak symmetry breaking (the focus point mechanism). In 
addition, as we will discuss in the next section, when m 3 / 2 ^ 4 TeV, the predicted relic 
density of the LSP well agrees with the currently observed dark matter density. 

Before closing this section, some discussion on the constraints from the gravitino prob- 
lem may be relevant. For the gravitino mass we are interested, overproduction of D gives 
the most stringent upper bound on Tr. Ref. |3] used averaged value of the observational 
values of D/H to set the bounds. If we adopt the largest value of observed D/H, upper 
bound on Tr becomes larger by factor 2 — 3. If so, the non-thermal leptogenesis may be 
consistent with the gravitino mass smaller than 4 TeV. However, notice that, in deriving 
Eq. (jl2j) . we considered the extreme case where Mi = |m$ and S c s = 1; in particular, if 
Mi is smaller, higher reheating temperature is required to generate large enough baryon 
asymmetry. 2 Thus, the non-thermal leptogenesis scenario is severely constrained from the 
point of view of the gravitino problem and the large gravitino mass is preferred to solve 
the conflict. 

3 Dark matter in LCSUGRA 

In the previous section, we have seen that the LCSUGRA scenario is good for the viable 
scenario of the non-thermal leptogenesis. In this section, we discuss that the LCSUGRA 
scenario can provide, in the reasonable portion of the parameter space, the dark-matter 
density consistent with the WMAP result. 

Before discussing the behavior of the dark-matter density, we first summarize the 
free parameters in our analysis. In the LCSUGRA scenario, all the higher- dimensional 
operators are assumed to be suppressed by the inverse powers of the cutoff scale M* 
much larger than the reduced Planck scale. Thus, as usual, all scalar particles obtain 
a universal SUSY breaking mass m through the effectively minimal Kahler potential, 3 

2 In addition, the primary purpose of Ref. |3] was to derive conservative constraints taking account 
of only the effects which are well understood. Thus, the hadrodissociations induced by the colored 
superparticles emitted from the gravitino decay were not included since such processes are very hard to 
estimate. Once they are included, the upper bound on Tr will become severer. 

3 There may be non-diagonal corrections to scalar mass squared from higher dimensional terms in 
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which is equal to the gravitino mass The universal gaugino mass min is obtained 

through higher dimensional operators in the superpotential which is suppressed by M* 
(where we assume unification of the SU(3) x SU{2) x £7(1) gauge group). The SUSY- 
breaking (scalar) 3 couplings, so called A parameters, are also generated through higher 
dimensional operators in the superpotential and they are expected to be of order the 
gaugino mass. 4 In the present analysis, we simply take A = 0, since their contributions 
to the following discussions are sub-dominant as long as \A\ <C m . In addition, we can 
obtain /i-parameter via the Giudice-Masiero mechanism [T7]. Then, the //-parameter is 
also suppressed by the large cutoff scale and becomes of the same order of the gaugino 
mass. Notice that the hierarchy /i <C m is consistent with the radiative electroweak 
symmetry breaking due to the focus point mechanism ;2J. Notably, in this case, the 
so-called 5-parameter is not a free parameter and it is equal to the gravitino mass (i.e., 
B = m ) for m ^> A. We consider that all of the above parameters are given at the GUT 
scale in the present analysis. Then, we have three parameters mo, mi/ 2 , [i (m ^> m^, n) 
at the GUT scale. 

At the electroweak breaking scale, the gaugino masses mi,m 2 and m 3 can be deter- 
mined by solving the RG equations, and are roughly given by mi ~ 0.4mi/ 2 , m 2 — 0.8m,!/ 2 
and m,3 ~ 2.8mi/2, respectively. Notice that the RG evolution of /i from the electroweak 
scale to the GUT scale is negligible, and hence, in the followings we take /i as the value 
at the electroweak breaking scale. We see that all the scalar mass parameters other than 
the one of the up-type Higgs scalar are much larger than gaugino masses and /i, thus, the 
gauginos and Higgsinos are much lighter than sfermions at the electroweak scale. Thus, 
the LSP is always lightest neutralino which is stable by i?-parity conservation. Such 
lightest neutralino becomes a good candidate of the dark matter as we see below. 

Postulating that the LSP becomes the dark matter, we calculate its relic density. Since 

heavier CP-even and CP-odd Higgs bosons as well as all the sfermions become much 

heavier than the neutralino, the leading processes in the neutralino pair annihilation are 

through s-channel exchanges of Z and CP-even Higgs bosons and also through the t- 

Kahler potential which are suppressed by the large cutoff M * . 

4 If the vacuum expectation value (VEV) of the field Z in the SUSY breaking sector is large as the 
reduced Planck scale Mg, we have the A terms of the order the gravitino mass to 3 / 2 - However, most of 
the dynamical SUSY breaking models suggest (Z) <C Mq |15M16| . Therefore, we consider |^4| ^ m 3 / 2 - 
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channel exchanges of the charginos and neutralinos |18j . In this case, tan/3 dependence 
of the relic density is small since the annihilation processes do not strongly depend on 
tan/3 except for the CP-even Higgs exchanges. 5 On the other hand, since we consider the 
situation \fj,\ ~ wii/2, the lightest neutralino may have a sizable fraction of the Higgsino 
components [H] and the relic density is sensitive to /i and mi/ 2 through the Higgsino 
fraction of the lightest neutralino. 

In Fig. H we show the parameter region in the (/i, m 2 ) plane where the neutralino 
density is consistent with WMAP result [131, ^dm/i 2 = 0.1126io:oi8i- 

In this computation, 

we have used micrOMEGAs 1.3.0 code [20J which includes all the possible co-annihilation 
effects. In the figure, the red points are for tan (3 = 5 and the blue ones for tan (3 = 30. As 
we mentioned, the result is insensitive to tan f3, and we find that the relic density is mainly 
determined by the parameters fi and rri2. It should be noted that the result in (/i, m^) 
plane is free from uncertainties arising from the analysis in the focus point region which 
we will discuss later. Especially, for ni2 ^ 350 GeV, tan (3 dependence of the relic density 
is negligible, since the neutralino becomes heavy enough so that the annihilation process 
is dominated by the channel into top-pair via Z-boson exchange. 6 The WMAP value 
of the relic density is realized on the line ~ 0.6 m2, and by comparing the bino mass 
parameter mi ~ 0.5 ni2, we find that the neutralino has a significant Higgsino fraction 
which makes annihilation more efficient. 

We next reinterpret the above results by using the LCSUGRA parameters. In the 
LCSUGRA, mo and tan (3 are uniquely determined for a given value of (/i,m 2 ), since the 
SUSY-breaking Higgs mixing parameter B is determined as B = m at the GUT scale. 

In order to understand how tan/3 depends on the LCSUGRA parameters, it is in- 
structive to see the tree-level condition for the electroweak symmetry breaking although 
important radiative corrections to the Higgs potential are taken into account in our numer- 
ical calculations. The tree-level minimization conditions of the effective Higgs potential 

5 tan/3 = v u /vd is the ratio of the two VEV's of the neutral Higgs fields, v Ul d = (^H° 

6 The processes into other lighter fermion pairs are much suppressed, since the amplitudes are propor- 
tional to the fermion masses which are required for helicity flip in s-wave processes. 
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(b) m top = 178 GeV 
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Figure 1: Dark matter constraints after WMAP (Q DM h 2 = 0.1126±g;gl|i) in the (£t,m 2 ) 
plane for (a) m top = 174GeV and (b) mt op = 178GeV. The reds and blue points correspond 
to tan/3 = 5 and tan/3 = 30, respectively. The black shaded regions are excluded by the 
chargino mass limit, m x ± > 104GeV 



are given by 



-rrir 



m 2 H tan 2 /? 



tan 2 p - 1 



sin 2/3 2 2 . oi ,2n 
— — {m Hu + m Hd +2\n\ ) 



(13) 
(14) 



where m 2 H denote the soft masses squared of up-type and down-type Higgs bosons, and 
all the parameters are evaluated at the electroweak scale. 7 Since all the parameters in 
Eqs. (jTH|) and (|T3j) are determined from three input parameters (m , mi/2, and //o, since 
B = mo), we can determine (mo, tan (3) for a given (/i,m 2 ) by solving above conditions. 

In order to obtain explicit relations between (mo, tan /3) and (/i, 777,2), we express the 
parameters in Eqs. (|13|) and (|14j) in terms of the input parameters as 



777 



777 



H u 
2 

B 



a u ml + b u |777 1/2 | 2 , 
a d ml + b d |777 1/2 | 2 , 
B = m , 



(15) 
(16) 
(17) 



7 In the actual numerical calculation, we use the running parameters at the typical stop mass scale, 
where the one-loop corrections to the Higgs potential tend to be small |22| . 
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where the coefficients ai and hi (i = u, d) are scale-dependent functions of dimensionless 
gauge and Yukawa coupling constants. Here, B ~ B is the consequence of the hierar- 
chical spectrum (mo 3> m\/2, A) of the LCSUGRA scenario. As discussed in Refs. 0, at 
the electroweak scale, a u is of order — 10 _1 while ~ 1 for m top ~ 170 — 180 GeV, 
which comes from universality of the scalar masses. The parameter bd becomes positive 
(0(0.1)) at the electroweak scale, which mainly comes from positive wino and bino contri- 
butions. On the other hand, b u becomes negative (0(1)) at the electroweak scale because 
of the positive gluino contributions to the stop mass squared which affects m 2 H through 
top Yukawa interaction. Substituting Eqs. (fT3j) . ()16|) and (fTTj) into the minimization 
conditions, we obtain explicit relation between (m ,tan/?) and ( / u,m 2 ) as 



These relations show that, when a u ~ 10~ 2 , the LCSUGRA scenario is consistent with 
the electroweak symmetry breaking for tan/3^ 10 [1. By using Eqs. (|18|) and (|19|). we 
obtain the WMAP constraint in the (tan (3, tuq) plane from the one in the (n, m 2 ) plane. 

With a more detailed numerical calculations, we relate the input parameters to the 
parameters at the electroweak scale, and derive the WMAP constraint in the (tan (3, m ) 
plane. The results are shown in Fig. which is converted from the one in Fig. In 
our numerical calculations, we have used the IS A JET 7.69 code [23] which takes into 
account the one-loop corrections to the effective Higgs potential and the two-loop RG 
evolutions of parameters. 9 In the figure, there is an upper bound on tan f3, which is 
from the lower bounds on the and m.2 in Fig. ^ The upper bound corresponds to 
parameters (tan/5, m , //, m 2 ) ^ (15,2 TeV, 140 GeV, 160 GeV) for m top = 174 GeV, 
and (25,4 TeV, 140 GeV, 160 GeV) for m top = 178 GeV. (The lower bound on the tan/3 
comes from the upper bound on m,2 1 TeV where we confine our attention.) From the 
figure, we find that the relic density is consistent for the WMAP result for mo ^ 2 TeV 

8 In Fig. |21 we reassure that the heavier CP-even, CP-odd Higgs bosons and all sfermions are much 
heavier than the neutralino, which is important to obtain the results in Fig. ^ 

9 It should be noted that, the lines in the figures show rough fitting of the results, since the code 
becomes somewhat unstable for mo 3> mi/2- 



tan/3 



m , 



(18) 




(19) 
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Figure 2: Cosmologically allowed regions of the relic density in the (tan (3, m ) plane for 
(a) m top = 174GeV and for (b) m top = 178GeV. 

for m top = 174 GeV and m ^ 4 TeV for m top = 178 GeV. 

As we have seen in the previous section, we find that the non-thermal leptogenesis and 
the gravitino problem require m 3 / 2 ^ 4 TeV. Using the fact that m = m 3 / 2 in LCSUGRA, 
we can see that the LCSUGRA scenario with a relatively large gravitino mass provides 
a consistent cosmological scenario which explains both baryon and dark matter densities 
in the present universe. 

Finally, we comment on the uncertainties in our results which comes from the technical 
difficulties in the precise calculation of the superparticle mass spectrum in LCSUGRA. 
Since some of the superparticles become much heavier than the electroweak-symmetry 
breaking scale in the focus point parameter region, an accurate determination of the mass 
spectrum requires careful treatments of various corrections. As a result, the numerical 
results are somewhat different from code to code |2*H I2*H] and our results may be affected 
by such uncertainty. To see its effect, we have repeated our analysis with SOFTSUSY 1.9 
code [2E]; we have find that the changes of mo and tan/? for a given set of (fi, 777-2) are 
at most ~ 30 %. Thus, we believe that our main conclusion, suggesting the LCSUGRA 
scenario with a relatively large gravitino mass, does not change although our quantitative 
estimations may be slightly changed. 
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4 Discussion 



In this letter, we have proposed a consistent cosmological hypothesis based on the LC- 
SUGRA scenario of the supersymmetry breaking. In the LCSUGRA hypothesis, the 
gravitino (as well as the sfermions) becomes as heavy as a few TeV without conflicting 
the naturalness of the electroweak symmetry breaking. If so, baryon asymmetry of the 
universe can be explained by the non-thermal leptogenesis via an inflaton decay without 
spoiling the success of the big-bang nucleosynthesis. We have also shown that, in the pa- 
rameter region suggested from the successful non-thermal leptogenesis, the relic density 
of the LSP explains naturally the dark matter density required from the WMAP results. 
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